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^! Abstract 

We present a detailed analysis on non-standard neutrino interactions (NSI) with electrons in- 
cluding all muon and electron (anti)-neutrino data from existing accelerators and reactors, in 
conjunction with the "neutrino counting" data (e^e~ -^ uv^) from the four LEP collaborations. 
i-£^ ' First we perform a one-parameter-at-a-time analysis, showing how most constraints improve with 

respect to previous results reported in the literature. We also present more robust results where the 



NSI parameters are allowed to vary freely in the analysis. We show the importance of combining 
LEP data with the other experiments in removing degeneracies in the global analysis constraining 
flavor-conserving NSI parameters which, at 90% and 95% C.L., must lie within unique allowed 
regions. Despite such improved constraints, there is still substantial room for improvement, posing 
a big challenge for upcoming experiments. 
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I. INTRODUCTION 



The historic discovery of neutrino oscillations constitutes the first evidence for physics 
beyond the Standard Model, and one would like to know to which direction it is pointing. 
Despite apretty good knowledge of the neutrino oscillation mechanism and the parameters 
involved pj] , the nature of the mechanism generatingmasses and mixings remains as elusive 
as ever. Neutrino mass models fall in various classes (21 involving models where neutrinos get 
„«=. a la s^a. aQaaflafl.Jio^wha.e neutnnos acquire „.... .adiat.valy due to the 



presence of extra Higgs bosons 



10 



JJ. 



and hybrid models, like those based on low en ergy 



14 



15|, 



16 



id 



supersymmetry with spontaneous (or bilinear) breaking of R-parity 
Starting with the seesaw 7| all those mechanisms carry with them modifications to the 
structure of the standard electroweak neutral and charged currents. The simplest of such 
modifications in the low energy regime may be written in the usual V — A form, similar to 
the four-Fermi interaction, characterized by some new coupling, which we generally call, in 
what follows, Non-Standard Interactions (NSI). Such interactions can arise in a variety of 
ways, e. g., from the exchange of Higgs and/or supersymmetric scalar bosons as well as a 
new heavy gauge boson, such as a Z'. Non-standard interactions can conserve flavor, their 
differences characterizing the violation of weak universality, known as Non- Universal (NU) 
NSI. Alternatively, they may violate flavor conservation, known as Flavor Changing (FC) 
NSI. 



On the other hand, current neutrino oscillation data, as inferred from the solar 18|, ll9|, l20 1 



and atmospheric 2l|, |22|, |23| neutrino experiments, leave significant room for the existence 



of sub-leading effects induced by NSI. In fact, NSI effects may be comparable to oscillation 
effects in solar neutrino physics, where indeed a new degenerate solution is still allowed 18|. 
At this stage laboratory experiments both from accelerators and reactors can play a crucial 
role, since the strongest sensitivity to NSI comes from this kind of experiments [2J]. Preci- 
sion measurements of oscillation parameters in long baseline oscillation experiments such as 
neutrino factories will also benefit from improved NSI studies. These could help resolvin, 
the confusion between the two found in Refs. 25 



26l | and further discussed in Refs. [271, |28| 



Here we address the current sensitivity on non-standard interactions as inferred from 



a global analysis of processes inyo 
have been performed in Refs. 



24 



vmg 



anti)-neutrinos and electrons. Previous analyses 



29|, Isd, l3l|- Our current analysis combines the relevant 
experimental "neutrino counting" data from e~^ + e^ —>■ u + u + '-/ obtained by the four LEP 
collaborations 

the t'e + e 



32 



33 



34 



35 



36 



37 



38 



39 



40 



4l| . and summarized in Ref. [42], with all 
i/g + e data obtained by LSND [43|], and the z/g + e — > z/g + e interaction studied 



in reactor experiments, namely: Irvine 4J], MUNU 45|| and Rovno 46|. For the muon 



neutrino case the relevant reactions are z/^ + e ^ z/^ + e and Ufj, + e -^ u^ + e, measured at 
CHARM II [471]. Our analysis is also novel in the sense that we adopt a model independent 
approach, as general as possible, allowing simultaneous variations of all NSI parameters. In 
particular, we compare the restrictions obtained varying only one parameter at a time, with 
those obtained in the case where all six flavor-conserving parameters are left free. 

The analysis sketched above is organized as follows: in Sec. [IT] the NSI are introduced 
and the relevant cross sections are expressed as function of NSI parameters, in Sees. Illll and 
IIVI we briefly present the data and the details of the x^ analysis. The results are presented 
in Sec. |V] and more discussion and outlook are given in Sec. IVD 



II. NON-STANDARD INTERACTIONS AND RELEVANT CROSS SECTIONS 



Neutrino NSI constitute an unavoidable characteristic feature of gauge models of neutrino 
mass, for example those where they arise from the admixture of isodoublet and isosinglet 
neutral leptons, like models of the generic seesaw type |7| . Typically the masses of the light 
neutrinos are obtained by diagonalizing the mass matrix 



Ml D 
D^ M 



R 



(1) 



in the basis z/, z/'^, where D is the standard SU{2) f/(l) breaking Dirac mass term, and 
Mr = M^ is the large isosinglet Majorana mass. In the absence of the isotriplet (type-II) 
Mluu term [f\ the scheme is called type-I seesaw 3|, 14, la, |6|. In models with spontaneous 
breaking of lepton number symmetry one has M^ oc 1/Mr, a feature that comes from the 
study of the scalar potential and holds both in the case of left-right models (gauged lepton 
number) [9i] and the case of majoron models (ungauged lepton number) ISj]. 

The structure of the associated effective SU{2) ® U{1) weak V — A currents is rather 
complex and deviates from standard ia\ . The flrst point to notice is that the heavy isosinglets 
will mix with the ordinary isodoublet neutrinos in the charged current weak interaction. As 
a result, the mixing matrix describing the charged leptonic weak interaction is a rectangular 
matrix K [71] which may be decomposed as 



K={Kl,Kh] 



(2) 



where Ki and Kh are 3x3 matrices. Note that the "effective" lepton mixing matrix K^ 



relevant in oscillation studies is non-unitary 
of such non-unitary effects see for example 



SOU . For papers addressing possible future tests 
5l| and references therein. The corresponding 



neutral weak interactions are described by a non-trivial matrix [3] K^K 



C 



w 



-Z^vl'^^K^Kvl . 



(3) 



2 sin 9w 

Such structure of the charged and neutral weak currents provides an origin for neutrino 
NSI. Note, however, that the smallness of neutrino mass, which follows due to the seesaw 
mechanism M^eS = M^ — DM^^D^ and the condition Ml <C Mr, implies that, barring 
fine-tuning, the magnitude of neutrino NSI and its effects are expected to be negligible. 

However this need not be so in general. Since the number m of SU{2) ® U{1) singlets 
is arbitrary, one may, for example, extend the lepton sector of the SU{2) (g) f/(l) theory by 
adding a set of two 2-component isosinglet neutral fermions, denoted u'^i and Si, in each 
generation. In such m = 6 models one can consider the 9x9 mass matrix 48|, |49(] 






D 





D^ 





M 





M^ 


/^ 



(4) 



(in the basis z/, u'^, S). The Majorana masses for the neutrinos are determined from 

Ml = DM-^fiM^'^D^ . 



(5) 



Since in the limit /i ^ the exact lepton number symmetry is recovered and neutrinos 
become massless 48| this scheme is sometimes called "inverse seesaw" 661 ]. 

This provides an elegant way to generate neutrino masses without a super-heavy scale, 
the smallness of the neutrino mass indicated by the oscillation interpretation of solar and 
atmospheric neutrino data is ascribed to the smallness of fi, which is natural in 't Hofft's 
sense: the symmetry of the theory is enhanced in the limit of vanishing /i. This automatically 
allows for a sizeable magnitude of neutrino NSI strengths, unconstrained by the smallness 
of neutrino masses [67| . 

The NSI which are engendered in this case will necessarily affect neutrino propagation 



52 



53| . They may 



properties in matter, an effect that may be resonant in certain cases 50, 
also be large enough as to produce effects in the laboratory. 

An alternative way to induce neutrino NSI is in the context of low-energy supersym- 
metry without R-parity conservation 



54 



55 



56| where one may also have, in addition to 



bihn^a. HQE Q also tn.i,... isolating couplings ,. the super-potential such as 



■^ijk-l^i-l^j-t^k 



(6) 
(7) 



where L, Q, E"^ and V^ are (chiral) super-fields which contain the usual lepton and quark 
SU{2) doublets and singlets, respectively, and i,j, k are generation indices. The couplings in 
Eq. (jH]) give rise at low energy to the following four-fermion effective Lagrangian for neutrino 
interactions with d-quark including 

Ces = -2V2Gf ^ ^af3 '^Lal'^^Lp dRl^d^ a,P = e,fx,T, (8) 

where the parameters ^ai3 represent the strength of the effective interactions normalized to 
the Fermi constant Gp- One can identify explicitly, for example, the following non-standard 
flavor- conserving NSl couplings 

IX' |2 

£ = V ' '^'^' (9) 



and the FC coupling 






where m^^.^ are the masses of the exchanged squarks and j = 1,2,3 denotes di,SL,bL, 
respectively. The existence of effective neutral current interactions contributing to the neu- 
trino scattering off (i-quarks in matter, provides new flavor-conserving as well as flavor- 
changing terms for the matter potentials of neutrinos. Such NSI are directly relevant for 



22 



23|. 



solar 18|, |l9|, |20|] and atmospheric neutrino propagation 21 , 

In what follows we consider a more general class of non-standard interactions described 
via the effective four fermion Lagrangian, 

- ^N?i = eipV2Gp{iy^^,LuMYPf) , (12) 

where Gp is the Fermi constant and e^^ parametrize the strength of the NSI. For laboratory 
experiments / is a first generation SM fermion {e,u or d). Here we analyze only processes 
involving electrons, so that in what follows we have only / = e. The chiral projectors P 
denote {L, R = {1 ± 7^)/2}, while a and (3 denote the three neutrino flavors: e, /i and r. 

In total there are 12 relevant parameters given by e^^. In order to constrain these we use 
experimental data reported by LEP (e+e" —>■ z/z/7), LSND (z/gC —>■ z/gc), CHARM II (z/^e or 
z/^e scattering) and reactor experiments (z/ge -^ Vf.e). The cross sections for the interactions 
of each experiment are given next. 



A. LEP cross section 



The e^e -^ uwy cross section can be calculated at tree level using the 'radiator' approx- 
imation to describe the photon emission 57| as 



a 



theo / 
LEP I 



dx I dc^ H(x,s^\s) ao^°°(s) , 



^, J-y, 



(13) 



where s is the center of mass energy, x = 2£'^/v^, E^ is the photon energy, o"o'^'^° = ctq^+ct^^^ 



is the cross section for the process e'^e 
H is given by 



uu and s = (1 — x)s. Here the 'radiator' function 



H{x, s^, s) 



2a 



TCXS^ 



1- 



X 



+ 



X^C^"! 



(14) 



where s^ = sm9-y and c^ = 1 — s^, with 6'^ being the photon emission angle. 

Working in the limit of vanishing W — j interactions but considering finite distance effects 
for the W propagator, the Standard Model e+e^ — > uu cross section is given as 



0-, 



SM, 



s) 



Gtt 



Mligl + gl) 



[{s - Mlf + {MzTz) 



w 



s + 2Ml 



w 



M^ fs + Ml 



2s 



W\^^^is + M^ 



9l- 



Ml{s - Ml) 



s - M|)2 + [MzTzY 



[s + M^^Y 



log 



Ml 

s + Ml 

Ml 



Ml 



(15) 



where N^, is the number of neutrino families, qr and ql are the SM electron coupling con- 
stants to the Z-boson, Mw, Mz and Tz are the W and Z-boson masses, and total Z decay 
width respectively. 



rNSI 

'o 



cr, 



The NU and FC components of the nonstandard cross section, ctq — ^q 
given by [29|: 



NU 



FC 
'O ' 



^0^"(^) 



cr, 



FC/ 



y ^s 



Qtt 



1^2 



+ 



\^ rvrv J I \^ nrv J 



2(nr^ I , r^ ) ^^js - M|) 



IT 



ee^^-^W 



[S + M2 ^2 



•■Wj 



log 



s + Mf 



w 



Ml 



G 



E ^s[[el,f + [e%f]. 



a^f3=e,^i,T 



Qtt 



Ml 



(16) 
(17) 



B. LSND and reactors cross sections 

The differential cross section for z/gC scattering processes in tlie presence of NSI can be 

written as 

+ (^^+5Zl^-r)(l-£) -(Mif + J^kaelkfeD'^el^], (18) 

wliere T is tlie electron recoil energy, vrtf^ is the electron mass and Ey is the incident neutrino 
energy. The effective couplings (jr^l are given as g^ = gR + e^^ and Ql = Ql + £ee- For the 
case of reactors, we have to exchange L hj R and vice-versa. 



C. CHARM II cross section 

In the presence of NSI the differential u^e — >■ UaC cross section relevant for the case of the 
CHARM II experiment is given as 



A = ^-^ 

ay 7r 




(i-yf 



(19) 



where ql^r = qlr + e^^^ and y = {\ — cos0*)/2 is called the inelasticity. Here 6* is the 
center of mass scattering angle. For the anti-neutrino case, we simply have to exchange L 
by R and vice-versa. 

III. THE DATA 

In what follows we will mainly focus on the effect of the six flavor conserving non-standard 
interactions in the above processes. Generalizing to include also the six flavor changing NSI 
parameters is straightforward but technically more complex and somewhat less motivated in 
view of strong, albeit indirect, bounds that follow from searches for lepton flavor violation. 
Let us now first briefly describe the relevant data used in our global analysis. 

A. The LEP data 

Neutrino-electron NSI will contribute to the cross section of the interaction e^e^ -^ z/z/7 
by increasing or decreasing the expected number of events. The best data on such interaction 



has been collected by the four LEP experiments: OPAL, ALEPH, L3 and DELPHI [32, 



33l. I34J . I35I . ISd I37I . ISSl . l39l . kd . |4l|. The reported measurements are compiled in Table 
H] 42j. The center of mass energy and luminosity for each of the four LEP experiments are 
given in the second and third columns of Table [B The background subtracted experimental 
cross sections and the Monte Carlo expectations are given in picobarns in columns four 
and five respectively. Column six reports the number of events observed after background 
subtraction. The efficiency e is given in column seven and finally, the last two columns 



report the kinematical cuts: x = E^/Ebeam, xt = xsin6'^ with 6^ t 



le angle between the 



photon momentum and the beam direction, and y = cos^^. As in 42|, we have found that 
our calculation for the Standard Model LEP cross section, Eq. ( TT3l) without including the 
effects of the NSI, disagrees with the Monte Carlo results quoted by the four collaborations. 
This might be due to additional specific experimental cuts beside the ones quoted in the 
last two columns in Table [H Such disagreements are included as an additional theoretical 
uncertainty which we have added in quadrature in the calculation of our errors. 

In total the four LEP experiments lead to 25 observables. Because of the small systematic 
error they have, we can assume that all of them are independent with no correlation between 
them. 



B. The LSND and reactors data 



The best measurements of the cross section for the z/gC and z/gC scattering processes 
have been performed in terrestrial experiments. The cross section for the elastic scattering 
interaction z/g + e~ ^ u^ + e~ was measured by the Liquid Scintillator Neutrino Detector 
(LSND) using a yU"^ decay-at-rest z/g beam at the Los Alamos Neutron Science Center. The 
detector is an approximately cylindrical tank containing 167 tons of liquid scintillator and 
viewed by 1220 photomultiplier tubes. The final neutrino-electron cross section is reported 
in Table [IT] ^4^ 



The Irvine 4J], the most recent MUNU 45|] and the Rovno 46| experiments have mea- 
sured the z/gC scattering by using neutrinos from reactors. The measured cross section is also 
reported in Table [TTl We also quoted the number of events for each experiment in column 
three and the recoil electron energy range in column two. 



C. The CHARM II data 



The CHARM collaboration used a massive 692 ton target calorimeter followed by a 



muon spectrometer to detect the z/^ + e ^ z/^ + e and z^» + e 



e scattering processes. 



The neutrinos were produced by a 450 GeV proton beam accelerated in the Super Proton 



TABLE I: Summary of 
the W~^W~ production 
second is systematic. 



the ALEPH, DELPHI, L3 and OPAL experimental data, collected above 
threshold. Wherever a double error is listed, the first is statistical and the 





v^ (GcV) 


-C(pb-i) 


a'^'^'' (pb) 


<tMC (pb) 


Nobs 


e{%) 


E-y (GeV) 


\y\ 


ALEPH 


161 


11.1 


5.3±0.8±0.2 


5.81±0.03 


41 


70 


XT > 0.075 


<0.95 


[3^ 


172 


10.6 


4.7±0.8±0.2 


4.85±0.04 


36 


72 


XT > 0.075 


<0.95 


[3^ 


183 


58.5 


4.32±0.31±0.13 


4.15±0.03 


195 


77 


XT > 0.075 


<0.95 




189 


173.6 


3.43±0.16±0.06 


3.48±0.05 


484 










192 


28.9 


3.47±0.39±0.06 


3.23±0.05 


81 










196 


79.9 


3.03±0.22±0.06 


3.26±0.05 


197 








]M 


200 


87.0 


3.23±0.21±0.06 


3.12±0.05 


231 


81.5 


XT > 0.075 


<0.95 




202 


44.4 


2.99±0.29±0.05 


3.07±0.05 


110 










205 


79.5 


2.84±0.21±0.05 


2.93±0.05 


182 










207 


134.3 


2.67±0.16±0.05 


2.80±0.05 


292 








DELPHI [35] 


















HPC 


189 


154.7 


1.80±0.15±0.14 


1.97 


146 


51" 


X > 0.06 


< 0.70 


FEMC 


183 


49.2 


2.33±0.31±0.18 


2.08 


65 


54" 


a;>0.2 


>0.85 


FEMC 


189 


157.7 


1.89±0.16±0.15 


1.94 


155 


50" 


a;<0.9 


<0.98 


L3 


161 


10.7 


6.75±0.91±0.18 


6.26±0.12 


57 


80.5 


> 10 


<0.73 


\M 














and 






172 


10.2 


6.12±0.89±0.14 


5.61±0.10 


49 


80.7 


Et>6 


0.80-0.97 


[sa 


183 


55.3 


5.36±0.39±0.10 


5.62±0.10 


195 


65.4 


> 5 
and 


<0.73 


[3^ 


189 


176.4 


5.25±0.22±0.07 


5.29±0.06 


572 


60.8 


Et>5 


0.81-0.97 


OPAL 


130 


2.3 


10.0 ±2.3 ±0.4 


13.48 ±0.22' 


19 


81.6 


XT > 0.05 


<0.82 


m 














or 






136 


2.59 


16.3 ±2.8 ±0.7 


11.30 ±0.20'' 


34 


79.7 


XT > 0.1 


< 0.966 




130 


2.35 


11.6 ±2.5 ±0.4 


14.26 ± 0.06 *" 


21 


77.0 






m 














XT > 0.05 


< 0.966 




136 


3.37 


14.9 ±2.4 ±0.5 


11.95 ±0.07'' 


39 


77.5 








161 


9.89 


5.3±0.8±0.2 


6.49±0.08'' 


40 


75.2 


a::T>0.05 


< 0.82 


m 














or 






172 


10.28 


5.5±0.8±0.2 


5.53 ±0.08'' 


45 


77.9 


a::T>0.1 


< 0.966 


m 


183 


54.5 


4.71± 0.34±0.16 


4.98±0.02'' 


191 


74.2 


a;T>0.05 


< 0.966 


[4JJ 


189 


177.3 


4.35±0.17±0.09 


4.66±0.03 


643 


82.1 


a;T>0.05 


< 0.966 



"Estimated from the Monte Carlo cross sections and the expected numbers of events. 
''Calculated from the expected number of events as predicted by the KORALZ event generator. 



TABLE II: Experimental measurements of the V(,e and V(,e scattering cross sections 



Experiment Energy range (MeV) Events 



Measurement 



LSND Ve 


- e 


10-50 


Irvine Ve - 


- e 


1.5- 3.0 


Irvine Ve - 


- e 


3.0- 4.5 


Rovno i>e 


— e 


0.6 - 2.0 


MUNU Ve 


— e 


0.7- 2.0 



191 a = [10.1 ± 1.5] X E^^MeM) x 10"^^ cm^ 

381 a = [0.86 ± 0.25] x ay-A 

77 CT = [1.7 ±0.44] X o-y-A 

41 0- = (1.26 ± 0.62) X 10-^'^cmVfission 

68 1.07 ± 0.34 events day "^ 



Synchrotron (SPS) for 2.5x10^^ protons on target. Approximately 10*^ neutrino interactions 
occurred in the detector. Data collected from 1987-1991 were 2677 ±82 events for reaction 
^n + e -^ ^'^±e and 2752 ±88 events in the u beam [47|. There was a neutrino contamination, 
of approximately 10% of the flux, in the muon-antineutrino electron scattering. 

The CHARM collaboration used these data to determine the values of the SM qa and 
gv coupling constants. Because of the quadratic dependence on the coupling constants in 
the cross section formula given in Eq. (iT9l) there is a well known fourfold ambiguity in the 
determination of gv and gA- A similar ambiguity in determining gy and g"^ has been removed 



m 



471 ] by combining the u^e and u^e scattering data obtained by the CHARM detector with 



the forward-backward asymmetry (Aps) in the e'^e^ -^ e~^e~ scattering at LEP 58|]. Here 
we will obtain a similar result in the context of constraining neutrino NSIs but using the 
"neutrino counting" LEP data. 



IV. THE x^ ANALYSIS 

Once we have defined in the previous sections the cross sections for each of the processes 
under consideration, and we have introduced all the experimental measurements relevant 
for our analysis, we proceed to perform a x^ analysis. 

For the LEP data, we can obtain a theoretical estimate of the expected number of events 
jytheo £qj, each of the 25 observables by using Eq. flT3l) . The integration of the cross section 
has been performed with the experimental cuts reported in last two columns of Table [B We 
have used the reported luminosity (C) and efficiency (e) for each experiment. 

We define the corresponding x^ function as 

25 



Xlep 



E 



1=1 



(Nl 



theo 



A^^ 



obs\2 



(20) 



where N^^^ is reported in Tableland Aj is the corresponding error. In the SM limit, our 



10 



cross section computations agree within 8% with the LEP Monte Carlo results, except for 
L3, where we have found up to a 20% discrepancy. Therefore, we have allowed for an extra 
10% theoretical systematic error added in quadratures to all LEP experiments 68j. We have 
neglected all correlations since statistical and systematic errors are small. 
For the v^e and z/gC scattering processes, we define the y^ as 

/-^thco _ ^cxp\2 

i « 

where the a^^^ are given by the experimental measurements and Aj are the corresponding 
errors reported in Table [TTl while 0"*^*^° = af^ + of^^ are the theoretical expectations con- 
sidering the effects of NSI calculated via Eq. ( ITSl) . Details of the analysis for LSND and 
reactor experiments constraining NSI in i/gC and z/gC scattering have already been given in 



Ref. [31 1 . 



For CHARM II we calculated the number of events using the cross section from Eq. ( TT9l) . 
With the NSI parameters fixed to zero we defined a normalization constant to reproduce 
the number of events reported by the CHARM II collaboration. 

For CHARM II data we have used 

XCHARM = 2^ T2 (22) 

where one observable stands for i/^e scattering and the other one for i/^e. 
The global x^ is simply the sum of the individual ones, 

XtOT = XlEP + Xv^e + XcHARM • (23) 



V. RESULTS 

The cross section for e e~ -^ z/z/7 including NSI is sensitive to all twelve e^o parameters. 
On the one hand the scattering interactions v^e and v^e are sensitive to six parameters: ef^^, 
with a = e, /i, r. On the other hand the elastic scatterings v^e and z/^e are sensitive to the 
other six parameters: e^^- 

In order to obtain constraints on the relevant NSI parameters, we first follow the most 



popular approach adopted by the majority of authors [2J, |29|]. It consists on varying only 
one parameter at-a-time and fixing the remaining parameters to zero. This way we obtain 
bounds on the twelve NSI parameters. However such one-parameter-at-a-time analysis is 
fragile and might miss potential cancellations in the determination of the restrictions upon 
NSI strengths. 

11 



As a second step, we assume that the new physics induces mainly flavor-conserving, 
effective NU neutral current interactions, so that the only relevant parameters are the six 
e^'^, a = e, fi, t. This is reasonable in view of the relatively strong bounds on lepton flavor 
violating processes. 

A. One parameter at-a-time 

In this section we constrain the neutrino-electron NSI parameters varying them one pa- 
rameter at-a-time. Because each cross section is sensitive to different parameters, depending 
on the parameter under consideration, the number of total observables used in the analysis 
will change. Table UTTl shows the Xmin^ ^^^ number of degrees of freedom (d.o.f), the ratio 
Xmm/d.o.f, the allowed range for each of the twelve parameters obtained by our x^ analysis 
and the limits obtained by previous analyses. 

TABLE III: Constrains of neutrino-electron NSI parameters by varying only one parameter at-a- 
time. Improvements are obtained compared with previous analyses. Read text for details. 



d.o.f. 


2 
Arnm 


xLn/d.O.f. 


90% C.L. Allowed Region 


Previous Limit [24, 31] 


29 


24.8482 


0.8568 


-0.03 < e^e < 0.08 


-0.05 < eg < 0.1 {u, 


.e) 


29 


22.4742 


0.7750 


0.004 < eg < 0.151 


-0.04 < eg < 0.14 {v 


ee) 


26 


22.1308 


0.8512 


\ej;^\ < 0.03 


le^^l < 0.03 (u^e) 




26 


22.1315 


0.8512 


k^J<0.03 


\e^^\ < 0.03 {i^.e) 




24 


21.8927 


0.9122 


-0.46 < e^^ < 0.24 


-0.6 < e^^ < 0.4 (e+e- - 


-^ uuj) 


24 


21.9072 


0.9128 


-0.25 < e^^ < 0.43 


-0.4 < eg^ < 0.6 (e+e- - 


-> i/i"y) 


31 


22.8752 


0.7379 


\e^^\ < 0.13 


|4| ~ 5 X 10-4 (^ ^ 


3e) 


31 


24.9885 


0.8061 


141 < 0-13 


|4| c^ 5 X 10-4 (/i ^ 


3e) 


29 


22.3062 


0.7692 


kerl < 0.33 


\eg\ < 0.4 (uee) 




29 


22.2107 


0.7659 


0.05 < lef^l < 0.28 


\eg\ < 0.27 (i^ee) 




26 


22.1308 


0.8512 


k^.l<o.i 


|e^J < 0.1 (u^e) 




26 


22.1312 


0.8512 


|4| < 0.1 


|e^,| < 0.1 (u^e) 





One sees how the inclusion of the LEP data leads to an improvement in the constraints 
for most of the NU NSI parameters. For example, from the last column in Tab 



24 



e mil one can 

are now 



see how previous constraints on eg^ coming from LSND and reactor data 

superseded. Our analysis also improves previous constraints on e!^^^. The inclusion of LEP 

data also improves the limits for eg^. Note that a nonzero eg is favored in this analysis. 
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though this has no statistical significance. We can also see that the ratio Xmm/d-o.f. is close 
to unity in the majority of the cases, meaning that the x^ is a good statistical indicator. 

Note, however, that there is no improvement in the constraints for the parameters ej^'J^, 
since these are dominated by the CHARM II data and the restrictions from /x — i> 3e for the 



NU and FC non-standard neutrino interactions, respectively 24l |. 

Here a comment on FC NSI is in order. Clearly, if there are FC NSI on neutrinos one 
expects, by SU(2) gauge symmetry, that these will induce also FC on charged leptons, which 
are rather strongly constrained by the non-observation of the corresponding LFV processes 
such as /i ^ 67, /i ^ 3e, /i ^ e conversion in nuclei, r — > /lee, r -^ ftp, etc. Indeed, given 
the existence of the effective NSI operators one obtains, by "dressing" with weak gauge- 
boson exchange, a corresponding effective NSI operator involving only charged leptons. 
However the loop diverges logarithmically. In this case a precise prescription must be given 
in order to estimate the corrections since the effective interactions are nonrenormalizable, and 
therefore there will be a dependence on the cuttoff scale A at which the theory is supposed 



to be renormalizable [24]. While the corresponding logarithmic terms can be rigurously 
computed when the physics producing the NSI lies at a large scale, it is certainly not so 
when it lies at a relatively low scale. The latter is precisely the case which is most relevant 
phenomenologically, for example, schemes like the extended seesaw, broken R-parity or 
radiative models of neutrino mass [id, lUl, |l2| . In these cases there is no model-independent 



way to rigurously compute the magnitude of the induced NSI among charged leptons in 
terms of that among neutrinos. 

It follows that so far NSI involving neutrinos are not strongly constrained, hence the 
importance of the constraints reported in Table IIIII in contrast with LFV constraints these 
are robust. 



Before clo sing this section let us mention that constraints coming from solar 18|, Il9| . 



atmospheric [2l|, |2^, |23|, and MINOS [59j data can not be directly compared with the 
bounds obtained here since those do not probe directly the NSI parameters but only a 
combination of them which effectively affects neutrino propagation in matter. For example, 
for the solar case, the relevant quantities, e and e' are two effective parameters which, for 
^afi ~ O5 ^^6 related with the vectorial couplings by: 

£ = - sin 623 el e' = sin^ 623 e^r - ^l , (24) 

with e^^ = Sap^^ap- Moreover, instead of just the NSI with electrons, one should in general 
take into account also the possible non-standard interactions with u and d -type quarks so 
that the effective NSI parameter becomes: e^^ = Yl, f=u d e^ Lp^ f / ^e (with nf the density of 
fermions in the medium). This leaves substantial freedom to new NSI-induced effects. 
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B. Flavor-conserving non-universal NSI 

Barring a theory of flavor, there is no guidance on the structure of the effective four- Fermi 
weak interaction. Generically new physics will lead to the violation of universality as well as 
the violation of leptonic flavor. In view of the relatively strong constraints on lepton flavor 
violating processes it is reasonable, as already mentioned, to first consider the case of purely 
flavor- conserving non-standard interactions, in general non-universal. In this case the only 
relevant parameters for our analysis are the six NU e^'J^. 

LSND and neutrino reactor data have been used previously in order to constrain the NU 



NSI parameters 3l|- It was noted that due to the nature of the elastic neutrino-electron 
scattering there is a fourfold ambiguity in the determination of the NSI parameters. The 
same happens when the analysis is performed for the non-universal parameters entering 
the Uf^e and u^^e scattering in the CHARM experiment. This fourfold ambiguity is clearly 
seen in the first two panels in Fig. [H The colored regions in the first panel show the two- 
dimensional projections in the £ee~^ee plane arising from the LSND and reactor data, while 
the corresponding restriction on the relevant parameters e^'^ arising from the CHARM 
experiment is displayed in the second panel of Fig. [H Finally, the third panel shows the 
projection of the constraints coming from the LEP data only on the parameters £:^|^. The 
dashed ellipses in the panels indicate the projections of the constraints following from LEP 
data only. 

These plots clearly indicate the complementarity between the "inclusive" LEP e~^e~ — »• 
uu'j data with those from reactors & LSND in constraining electron-type NSIs. Similar 
complementarity holds between LEP data and CHARM II data when constraining muon- 
type NSIs. 

In the global analysis where Xtot is the addition of the LEP, CHARM, LSND and reactor 
pieces, one clearly sees how the above fourfold ambiguities are eliminated. This is illustrated 
in Fig. [2] where we show the allowed regions that arise from the global Xtot after taking a 
projection over two parameters. The shaded (colored) areas shows the 90 %, 95 %, and 99 
% C.L. allowed regions (corresponding to AxI-qt = ^xLin + 4.61, 5.99, 9.21 respectively). 
The constraints derived from this analysis are collected in Table IIVI and compared with 
the results discussed in the previous section. One can see that the interplay between the 
different experiments, namely, the combination of the LEP neutrino counting results with 
the remaining data, plays a crucial role in providing constraints almost as stringent as in 
those obtained in a one-at-a-time analysis. Needless to say the global analysis establishes 
the robustness of these constraints since we are allowing all the six parameters to vary. 
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FIG. 1: Flavor conserving NSI allowed by each experiment in our analysis. The colored (shaded) 
regions in the first and second panel show the allowed 90 %, 95 % and 99 % C.L. regions from the 
Ve& and CHARM II data, respectively. For the LEP data, we take the projection of the six free 
parameters, over the two displayed NSI parameters in each case. In this case, solid, dashed and 
dot-dashed lines show the 90 %, 95 % and 99 % CL allowed regions respectively. 

TABLE IV: Constraints for the flavor conserving parameters. We have included the data from 
reactor, CHARM II and LEP experiments and allowed the six flavor conserving NSI to be present. 
The global minimum for the \^ analysis is X^^in — 23.13, Xmin/ d.o.f. = 1.22. We show the 90 
% C.L. values obtained after taking a projection over two parameters. For comparison, we show 
the constraints for the case in which only one parameter is allowed to vary, and finally, we also 
compare with previous reported results for the case of one parameter at a time. 





90% C.L. Allowed Region 


One parameter 


Previous hmits 




-0.14 < ef-g < 0.09 


-0.03 < efg < 0.08 


-0.05 < efg < 0.1 


■-ee 


-0.03 < ef, < 0.18 


0.004 < ef, < 0.15 


0.04 < ef, < 0.14 




-0.033 < e^^ < 0.055 


\ejl^\ < 0.03 


le^^l < 0.03 




-0.040 < e^^ < 0.053 


le^^l < 0.03 


|e«^| < 0.03 


e^r 


-0.6 < e^^ < 0.4 


-0.5 < e^^ < 0.2 


krri < 0.5 


e?r 


-0.4 < e^^ < 0.6 


-0.3 < e^^ < 0.4 


|e?rl < 0.5 



VI. DISCUSSION AND CONCLUSIONS 

We have given a detailed analysis on non-standard neutrino interactions with electrons 
following from combining muon and electron (anti)-neutrino data collected in existing accel- 



15 



0.5 



I— 1 ~^ 



■0.5 



-1 1 1 1 1 1 1 1 r- 



-1 1 1 r- 




L,R 



_l I I I I I I I I I I I I I l_ 

-0.5 0.5 



R 

FIG. 2: Global analysis results for the flavor conserving NSI (details in Table lIVp . The plot 
shows the regions allowed at 90 %, 95 % and 99 % C.L. arising from the marginalization of the six 
parameters over different flavors e^ , for I = e,ji and r, as indicated. For each flavor the allowed 

r D 

region is unique, except for the case of e^'^ , where there is a second solution, which is allowed only 
at the 99% C.L. 

erators and reactors, with the high energy "neutrino counting" data from LEP. Except for 
e^'^ and most PC NSIs, the inclusion of the LEP data within a simple one-parameter-at-a- 
time analysis improves upon previous constraints on the flavor- conserving NSI parameters. 

Barring a fundamental theory of flavor, there is no theoretical guidance on the flavor 
structure of the NSI that presumably result from the basic underlying theory producing 
neutrino masses. As a result the expected modifications in muon and electron (anti)-neutrino 
interactions involve the various components of the NSIs. Given this, it is necessary to 
perform a more general and robust analysis in which ideally all NSI parameters are allowed 
to vary freely. As a first step we have considered the case of non-universal NSIs. Our results 
indicate a strong complementarity between the "neutrino counting" data and the rest in 
removing the ambiguous determination of NSI parameter bounds. We have obtained unique 
allowed regions at 90% and 95% C.L. in NSI parameter space. Our improved constraints 
still leave substantial room for improvement, posing a big challenge for the next generation 
of neutrino experiments. 
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